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Abstract 

We point out that the recent data from the Sudbury Neutrino Observatory, together with other 
relevant measurements from solar and reactor neutrino experiments, convincingly show that the 
flavor transitions of solar neutrinos are affected by Mikheyev-Smirnov-Wolfenstein (MSW) effects. 
More precisely, one can safely reject the null hypothesis of no MSW interaction energy in matter, 
despite the fact that the interaction amplitude (formally treated as a free parameter) is still weakly 
constrained by the current phenomenology. Such a constraint can be improved, however, by future 
data from the KamLAND experiment. In the standard MSW case, we also perform an updated 
analysis of two-family active oscillations of solar and reactor neutrinos. 

PACS numbers: 26.65.+t, 13.15.+g, 14.60.Pq, 95.55.Vj 
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I. INTRODUCTION 



The Sudbury Neutrino Observatory (SNO) experiment has recently released new data 
0] with enhanced sensitivity to neutral-current (NC) interactions of solar neutrinos in deu- 
terium. Charged current (CC) and elastic scattering (ES) events have also been statistically 
separated from NC events in a model- independent way, i.e., without using priors on the 
neutrino energy spectrum shape 0]. These data corroborate the explanation of the solar 
neutrino deficit 1^ in terms of (dominant) two-family z/g Ua flavor transitions = z/^,t), 
which have convincingly emerged from the combined data of previous solar neutrino exper- 
iments (Chlorine 0, Gallium ]1 H , Super-Kamiokande (SKi_0,l3, and SNO BHEil) 
and of long-baseline reactor oscillation searches at KamLAND Moreover, the new SNO 
data appear to forbid relatively high values of the neutrino mixing angle (close to maxi- 
mal mixing) and of the squared mass difference Sm? (close to the CHOOZ upper bound), 
which were marginally allowed prior to 5! (see, e.g., HIH). In the current global fit 
the mass-mixing parameters appear to be tightly confined in the so-called large mixing angle 
(LMA) region, and especially in a subregion previously denoted as LMA-I 

In the LMA parameter range, flavor transitions between Ve and should be significantly 
affected by the neutrino interaction energy difference V = Ve — Va arising in solar (and 
possibly Earth) background matter [l6l[l7|. 

V{x) = V2GfN,{x) , (1) 

where Ng is the electron number density at the point x. The associated flavor chan ge, known 
as Mikheyev-Smirnov-Wolfenstein (MSW) effect |l6|, should occur adiabatically flq in the 
solar matter, for LMA parameters. In the context of Hamiltonian (7i) evolution of 2u active 
flavors, the MSW effect enters through a dynamical term Tidyn in matter, in addition to the 
kinetic term Tikin in vacuum: 





(^dyn + ^kin) , (2) 



where 



and 




_Sm^ / - cos 2^12 sin 2^12 \ , . 

4E [ sin 2^12 cos 2^12 j ' 

E being the neutrino energy. 

In a previous work we pointed out that, while the evidence for Tikin 7^ was over- 
whelming, the phenomenological indications in favor of Tidjm 7^ (and thus of MSW effects) 
were not as compelling. In particular, we introduced in [19| a free parameter ausw modu- 
lating the overall amplitude of the dynamical term Tidyn through the substitution 

V ausw - V , (5) 

both in the Sun and in the Earth. We showed that ausw was poorl y co nstrained, despite an 
intriguing preference for the standard MSW expectation omsw ~ 1 • The null hypothesis 
flMsw = was not clearly disproved by any single experiment, and could be rejected at a 
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relevant confidence level (Ax^ ^ 13, formally equivalent to y/Ax^ ^ 3.5cr) only in the global 
fit. We concluded that the available phenomeno logical data clearly favored MSW effects in 
solar neutrinos, but did not prove unequivocally their occurrence. We deemed it necessary 
to wait for new KamLAND or SNO data, in order to clarify the situation and to probe MSW 
effects with higher statistical significance |l9l |. 

In this work, we point out that the recent SNO data ^ contribute significantly to dis- 
prove the null hypothesis of no MSW oscillations. In the global combination of solar and 
reactor data, we find that, with respect to the (preferred) standard case omsw ~ the 
null hypothesis cimsw = can be safely rejected at the level of ~ 5.6a, despite the fact the 
allowed range of omsw is still rather large. In other words, the evidence in favor of MSW 
effects is now very strong, although precision tests of the MSW physics cannot be performed 
until new, high statistics KamLAND data become available (as we show later). 

In the following sections, we analyze the current solar and reactor neutrino phenomenol- 
ogy with an increasing degree of dependence on assumptions about the MSW effect.^ In 
Sec. II we do not make any hypothesis about MSW effects, and show that SNO data alone, 
as well as a model-independent SNO+SK combination, constrain the energy-averaged 
survival probability (Pee) to be significantly smaller than 1/2. This fact, by itself, excludes 
the vacuum case amsw = (which would predict {Pee) > 1/2 in the LMA region selected by 
KamLAND), and proves that dynamical effects must occur in solar neutrino propagation 
with unspecified amplitude omsw > 0. In Sec. Ill we fit all the available solar and reac- 
tor data with (5m^, 9i2, omsw) taken as free parameters. We find that MSW effects with 
standard amplitude (omsw = 1) are favored, while the null hypothesis (cmsw = 0) can be 
safely rejected at the ~ 5.6a level. However, we show that the allowed range of omsw is still 
very large, and can be significantly narrowed only by future KamLAND data. Assuming 
standard MSW effects (omsw = 1), we perform in Sec. IV an updated analysis of the 2u 
kinematical parameters (5m^, sin^ O12). We briefly discuss the impact of mixing in Sec. V, 
and conclude our work in Sec. VI. 

II. MODEL-INDEPENDENT CONSTRAINTS 

It has been shown in fiO] (see also that the SK and SNO experiments probe the same 
energy-averaged Ue survival probability (Pee) to a good accuracy, provided that the detector 
thresholds are appropriately chosen. For the kinetic energy threshold (Tsno = 5.5 MeV) 
and energy resolution characterizing the latest SNO data l] , we find that the equivalent SK 
threshold is E^k — 7.8 MeV in total energy.^ For equalized thresholds, the SK ES flux and 
the SNO NC and CC fluxes are linked by the exact relations |2(| 



where r = 0.154 is the ratio of (properly averaged) u^^j. and z/g CC cross sections, and $b is 
the true flux from the Sun. From the above equations, one can (over) constrain both $5 



^ In any case, we assume active flavor oscillations only, and do not consider hypothetical sterile neutrinos. 
^ In practice, adopting this threshold amounts to discard the first three bins of the SK energy spectrum in 
the range E e [5, 8] MeV. 




<l>BKPee)+r(l-(Pee))] , 
^B{Pee) , 



(6) 
(7) 
(8) 
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and (-Pee) in a truly model-independent way, namely, without any prior assumption about 
the energy profile of Pee or about $b predictions in standard solar models (SSM). 

Figure 1 shows the current constraints on $b and on (Pee) as derived from the final SK 
ES data and from the latest SNO CC and NC fluxes [l[ (correlations included [l^)- The 
constraints are shown both by individual bands and by their combination at the 3a level 
(Ax^ = 9). The projections of the SNO+SK combination (solid ellipse in Fig. 1) provide 
the ranges 

$5 = (5.5 ± 1.2) X lO'' cm-^s"^ {3a) , (9) 
in good agreement with SSM predictions and 

(Pee) = 0.3ll°:J^ ^3^) _ (IQ) 

The above 3a limits on (Pee) are in very good agreement with the "3(t range" obtained by 
naively triplicating the errors of the SNO CC/NC flux ratio, which is a direct measurement 
of (Pee): ^^'^o/^^^o = 0.306 ± 0.105(3cr) Q. However, as emphasized in 0, the errors of 



the CC/NC ratio are not normally distributed, and should not be used in fits. Conversely, 
our bounds in Eq. ()10|) are statistically safe and well-defined, and will be used in the following 
discussion. 

The above SK+SNO constraints appear to be currently dominated by the SNO data, as 
shown by the dotted ellipse in Fig. 1. In particular, the upper bound on the z/g survival 
probability, 

(Pee) < 0.43 (3a) , (11) 

can be basically derived from the SNO (CC+NC) data |1|] alone. The upper limit in Eg. ()11|) 
is significantly stronger than the one derived in [ilj prior to the latest SNO data m. In 
particular, we have now robust, model-independent evidence that Pee is definitely smaller 
than 1/2 at > 3a level. This inequality has important consequences for both the dynamical 
and the kinematical term in Eq. (j2)). First, in the Sm"^ range accessible to KamLAND and 
below the CHOOZ bound (5m^ ~ 0(10"''^^) eV^), the absence of the dynamical MSW term 
'Hdyn (i-e., the case cmsw = would imply (Pee) > 1/2 (see, e.g, [3l), contrary to Eq. (fTT|) . 
Second, assuming standard MSW dynamics (amsw = 1)? the inequality in Eq. (fTTj) allows to 
place upper limits on the kinematical parameters and sin^ 6*12 (see, e.g., the discussions 
in til HQ). 

Figure 2 illustrates the previous statements through isolines of (Pee) in the (5m^, sin^ Ou) 
parameter range relevant to KamLAND, for both omsw = (left panel) and omsw = 1 
(right panel). The superposed grey region is allowed by the 3a model-independent bounds 
in Eq. (jTU)). No such region exists in the case omsw = 0, which is therefore rejected at the 
3a level (at least). In the standard MSW case (left panel), the allowed region appears to 
confine 6ni^ below ~ 2 x 10~^ eV^ and sin^6'i2 below ~ 0.4. In particular, the latest SNO 
data significantly contribute to reject maximal mixing (sin^^i2 = 1/2) and to reduce the 
likelihood of the so-called 3| LMA-II parameter region, as emphasized in 0]. 



Summarizing, the latest SNO CC and NC data [l|, either by themselves or in combination 
with the SK ES data j^, provide the strong, model-independent upper bound (Pee) < 0.43 at 
3a. In the context of 2^ mixing, and within the mass-mixing region probed by KamLAND, 
this bound allows to reject the null hypothesis (amsw = 0), and provides upper limits on 
the mass-mixing parameters in the standard MSW case (qmsw = !)• In the next Section, 
we examine the more general case of variable ausw, in order to test whether current and 
future data can significantly constrain, by themselves, the size of matter effects. 



4 



III. CONSTRAINTS ON THE MSW DYNAMICAL TERM 



In this section we present the results of a global analysis of solar and reactor (Kam- 
LAND+CHOOZ) data with ((5m^, sin^ ^12, answ) unconstrained. The latest SNO data III 
are incorporated according to the recommendations in i^l- The reader is referred to 
for other details of the analysis. 

Figure 3 shows the results of the global fit, in terms of the function Ax^(aMsw) after 
{Sm'^, sin^ 612) marginalization. Such marginalization is appropriate to test the size of Tidyn 
independently of Tikin- It can be seen that the best fit is intriguingly close to the standard 
case (ctMsw = 1), although there are other acceptable local minima over about three decades 
in aMsw- As discussed in [l9[ for the case of variable cimsw, the Sm^ range allowed by solar 
neutrino data sweeps through the tower of LMA-n solutions allowed by KamLAND, leading 
to a series of "bumps" in the Ax^ function (solid hne). Such features are unavoidable, as 
far as KamLAND allows multiple solutions in the mass-mixing parameter space. However, 
the situation should improve with higher KamLAND statistics. Assuming that KamLAND 
will confirm the current best-fit solution in the ((5m^, sin^ ^12, cimsw) space, and simulating 
the corresponding KamLAND data, we obtain the prospective dotted and dahed curves in 
Fig. 1, which refer, respectively, to a fivefold and tenfold increase of the present statistics 
(54 events ^3)- It appears that, with the help of a few hundreds KamLAND events, the 
global fit of solar and reactor data can pinpoint the predicted size of MSW effects within 
a factor of ~ 2, allowing future "precision tests" of this effects (e.g., to probe additional 
nonstandard interactions). 

Although the current bounds on omsw appear to be rather weak, the rejection of the null 
hypothesis omsw = is quite strong, and corresponds to a significance level of Ax^ — 32, i.e., 
~ 5.6(7. Summarizing the results of this and the previous section, we can state that current 
solar and reactor data reject the hypothesis of no MSW effect at > 5cr level, with a > 3a 
contribution from the recent SNO data . Therefore, in our opinion, the phenomenological 
indications in favor of MSW effects can now be promoted to the level of evidence. 



IV. CONSTRAINTS ON KINEMATICAL MASS-MIXING TERM 

In this section, assuming standard MSW dynamics, we update our previous bounds [U] 
on the mass- mixing parameters ((5m^, sin^ Ou) which govern the kinematical term Tikin- The 
reader is referred to [l^, 2ll| for technical details. Here we just add that the statistical 
correlations of recent SNO data are incorporated through a straightforward generalization 



of the pull approach [21[, as explicitly described in j26|. We have checked that our analysis 
"SNO data only" reproduces the results of [i[ with very good accuracy. Finally, we have 
updated the total rate and winter-summer asymmetry from Gallium experiments [27]. In 
total, we have 84 solar neutrino observables, plus 13 KamLAND bins. 

Figure 4 shows the results of our fit to all solar neutrino data.^ The best fit (Xmin = 72.9) 
is reached at = 5.7 x 10~^ and sin^ 9^ = 0.29. The upper and lower bounds on the 
mass-mixing parameters are in good agreement with the results in and confirm that the 



We used to add CHOOZ data |3| in order to strengthen the upper bound on 5m'^ However, 
current solar neutrino data make this addition no longer necessary in the context of mixing with 
standard MSW effects. 
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solar neutrino parameter space is steadily narrowing. 

Figure 5 incorporates the analysis of KaniLAND data fl2| as in The best fit (Xmin — 
79.7) is reached at 6m^ = 7.2 x 10"^ and sin^6'i2 = 0.29 (LMA-I solution), while the second 
best fit (LMA-II solution) is only marginally allowed at the Ax^ = 9.4 level (~ 99% CL. 
for A^DF = 2). Also in this case, we find good agreement with the results in modulo 
the obvious transformation from our linear abscissa sin^ 6*12 to their logarithmic abscissa 
tan^ 612. 

In conclusion, the kinematical 2i/ mass-mixing parameters appear to be strongly con- 
strained in a basically unique region (LMA-I), with only a marginal possibility left for the 
LMA-II region. The decrease of the previous LMA-II likelihood M is an important contri- 
bution of the latest SNO data 



V. COMMENTS ON THREE-FAMILY MIXING 

So far, we have assumed flavor oscillations in the active 2i/ channel Va {^a being a lin- 
ear combination of and Vr) driven by the {5m? ^ 612) parameters. The (z/^, Ur) combination 
orthogonal to z/q is probed by atmospheric z/^ — > oscillations, with different parameters 
(Am^,^23) ^281. As far as the third mixing angle is zero (and 5rn? / /S.rn? <C 1), the 
two oscillation channels are practically decoupled, and all our previous considerations hold 
without changes. However, for small but nonzero ^13, the 3z/ survival probability deviates 
from the 2v case for both solar and KamLAND z/g oscillations: 

P,3:^(l-2sin2^^i3)Pe':. (12) 

In particular, for cmsw = 0, the minimum value of {Pee) in the right panel of Fig. 2 can 
slightly decrease from 1/2 to 1/2 — sin^6'i3. Until very recently, the upper bound on 6*13 
(dominated by CHOOZ and atmospheric data) could be quoted as sin^ 6*13 < 0.05 (Scr) [IJ, 
2^, leading to Pe^e^('^MSw = 0) > 0.45. A new SK atmospheric data analysis js^l, however. 



appears to imply the weaker bound sin^^is < 0.067 (3cr) [31|, leading to -Pe^e'^(aMsw = 0) > 
0.43. In both cases, there is no overlap with the experimental upper bound in Eq. (Illj). 
Therefore, the null hypothesis omsw = can be rejected at the 3a level also in the 3z/ 
mixing case, using only SNO(+SK) data. 

In the more general case of variable cimsw; we have not performed the 3z/ generalization 
of the analysis in Sec. III. Our educated guess is that an allowance for small values of 61^ 
should only slightly weaken — but should not spoil — the main results discussed therein. 



VI. CONCLUSIONS 

In recent years, solar and reactor neutrino data have been shown to be consistent with 
(and to favor) Mikheyev-Smirnov-Wolfenstein effects in the flavor evolution of solar neutri- 
nos. However, in our opinion, the null hypothesis of "no MSW effect" could not be safely 
rejected [19]. The current situation appears, however, more favorable. In this work we 
have pointed out that recent SNO data strongly favor the occurrence of MSW effects 
in the solar matter and, together with world solar and reactor data, provide a many-sigma 
rejection of the null hypothesis. We have also performed an analysis where the MSW in- 
teraction energy is freely rescaled, and found waek constraints on the scaling parameter. 
These constraints can be potentially improved by higher-statistics KamLAND data, which 
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will then allow more precise tests of the MSW dynamics. In the standard MSW case, we 
have also performed an updated analysis of two-family active oscillations of solar and reactor 
neutrinos. 

We conclude by observing that, although MSW effects are an unavoidable consequence of 
the standard theory of electroweak interactions, their basic confirmation in the current neu- 
trino phenomenology represents an important and reassuring experimental accomplishment, 
which strengthen our confidence in the emerging picture of neutrino masses and mixings. 
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FIG. 1: Results of the model-independent analysis of SNO (CC and NC) and SK (ES) neutrino 
fluxes. The projections of the solid ellipse provide 3a bounds on the neutrino flux and on 
the energy-averaged survival probability (Pee)- 
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FIG. 2: Isolines of {Pee) with standard MSW effects (aMSW=i) and with no matter effect (answ = 
0). The gray region is allowed by the SK+SNO combination. No such region exist in the absence 
of MSW effects. 
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Bounds on o^sw (Solar + CHOOZ + KomLAND) 
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FIG. 3: Bounds on cmsw (considered as a continuous free parameter), including all current solar, 
CHOOZ, and KamLAND data (solid curve). Prospective KamLAND data with higher statistics 
are used to draw the dotted and dashed curves. See the text for details. 
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FIG. 4: Results of the 2j/ oscillation analysis of solar neutrino data, for standard MSW effects. 
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FIG. 5: Results of the 2i/ oscillation analysis of solar and KamLAND data, for standard MSW 
effects. 
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